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Abstract 

We present a phenomenological study of the rare double radiative decay 
B i^77 in the Standard Model (SM) and beyond. Using the operator prod- 
uct expansion (OPE) technique, we estimate the short-distance (SD) contri- 
bution to the decay amplitude in a region of the phase space which is around 
the point where all decay products have energy ~ m^/S in the rest frame of 
the i?-meson. At lowest order in 1/Q, where Q is of order uib, the B — > K'jj 
matrix element is then expressed in terms of the usual B ^ K form factors 
known from semileptonic rare decays. The integrated SD branching ratio in 
the SM in the OPE region turns out to be l^B{B Kj'y)^[/^ ~ 1 x IQ-^. 
We work out the di-photon invariant mass distribution with and without the 
resonant background through B — > K{r]c, Xco} ~^ Kj'j. In the SM, the reso- 
nance contribution is dominant in the region of phase space where the OPE 
is valid. The present experimental upper limit on Bg t~^t~ decays, which 
constrains the scalar /pseudoscalar Four- Fermi operators with t~^t~, leaves 
considerable room for new physics in the one-particle-irreducible contribution 
to S — > X77 decays. In this case, we find that the SD B — > K^^f branching 
ratio can be enhanced by one order of magnitude with respect to its SM value 
and the SD contribution can lie outside of the resonance peaks. 
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1 Introduction 



Rare S-decays mediated by flavor changing neutral currents (FCNC) provide a rich 
laboratory to study the flavor structure of the Standard Model (SM) and its exten- 
sions. Such processes are forbidden in the Born approximation and are generated 
by loops, which gives them a high sensitivity to new physics (NP). A prominent ex- 
ample is the both experimentally and theoretically well-studied 6 — > 57 transition, 
exemplifled by the inclusive B — > Xg^ and the exclusive B — > K*^ decays (for a 
detailed review, see [1]). 

Radiative 5-decays involving two photons, such as B ^ 77, B X'-f-f or 
77 decays, are of further interest. Among the Cabibbo-favored h ^ s 
transitions, the decays Bg — > 77 and B — > ^577 have been theoretically analyzed 
in the SM some time ago [2] and later on by taking into account the leading-order 
QCD corrections in the framework of the electroweak effective Hamiltonian [3-6]. 
Recalling that at the quark level the FCNC double-photon decay is mediated by 
h — > S77, it receives one-particle-reducible (IPR) contributions from the b ^ s'j 
transition plus an additional photon and a one-particle-irreducible (IPI) term from 
a fermion loop with the two photons emitted from that loop. The correlation with 
b — >■ S7 implies that NP, which predominantly induces contributions to the IPR 
piece, gives no drastic effects in the double radiative decays because of the strong 
constraints from B Xs'~f data. This has been shown, for example, for the two- 
Higgs-doublet model [7]. On the other hand, NP in the IPI contribution, that is 
in Four-Fermi operators, is subleasing in 6 — > 57 decays and hence can be sizable 
in 6 — s> S77 transitions where it enters at the same order as the IPR contributions. 
So far only experimental upper limits on double radiative decays have been set (at 
90% C.L.): 

B(Bd -f-f) < 1.7 X 10-^ (BaBar) [8] , 

B{Bg ^ 77) < 1.48 X 10-^ (L3) [9], (1) 

which are about two orders of magnitude above their respective SM values, e.g. [6]. 

While many theoretical investigations have been dedicated to Bg — > 77 and 
B — > Xs77, less attention has been devoted to cxchisivc B K^*'>'y'y decays. For 
instance, only two groups [10,11] have estimated the B — > ^^77 branching ratio, 
with the result given as ~ (0.5 — 5.6) x 10~^, depending on cuts on the photon 
energies. Note that B K77 decays have a complicated matrix element due to 
the non-local IPR contribution. In the aforementioned studies the IPR piece has 
been solely parametrized in terms of a long-distance contribution using hadronic 
models, i.e. vector meson dominance. There is also a sizable resonant background 
from charmonia through B — > K{cc) — > i^77. Since B — > ^^77 decays provide a 
reasonable experimental signal for not too soft photons and with branching ratios in 
the 10~^ range, they are in reach of the currently operating S-factories at SLAC and 
KEK and at experiments at the hadron colliders, Tcvatron and LHC. Therefore, a 
detailed study of these decays in the SM and beyond is of great interest, in particular 
due to the complementarity with 6 — > S7 decays with respect to NP searches. 
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Our main objectives in this paper are i to write down a model-independent de- 
scription of exclusive B Kjj decays and ii to find out to what extend one can 
extract short-distance (SD) physics from these decays, that is, how big is the SD 
contribution and can it be experimentally isolated from the resonance contribu- 
tions. The idea here is to use the freedom of a 3-body decay and choose the photon 
energies such that the virtualities of the intermediate quarks in the IPR b 577 di- 
agrams are hard. Then, these scales of order rri}, are integrated out with the effective 
field theory methods using heavy-quark-effective-theory (HQET) [12] and the soft- 
coUinear-effective-theory (SCET) [13,14]. The resulting operator product expansion 
(OPE) enables us to express the matrix element in terms of the B K form fac- 
tors in some region of the phase space. There are other mechanisms contributing 
to S — > ii'77 decays, such as when one photon gets emitted from the spectator 
quark or annihilation diagrams, which we comment on. We further incorporate the 
rjci Xco-resonance contributions using a phenomenological parametrization a la Breit- 
Wigner to see whether they can be distinguished from the SD ones obtained from 
the OPE. We work out possible NP effects in the B K'j'j spectra and branching 
ratios. In particular, we illustrate the impact of a scenario with enlarged (sfe)(r^r~) 
couplings, which substantially modifies the IPI b s'-ff contribution but does not 
violate existing constraints on any other rare 5-decays. 

The remainder of this paper is organized as follows: After a brief description of 
the effective Hamiltonian formalism for double radiative b s'-fj decays in Section 
12 we evaluate in Section El the (non-local) matrix element of -B — ii'77 decays with 
the OPE and present the resulting decay amplitude and distributions. Additional 
contributions to the B K'-yj amplitude, namely photon radiation off the spec- 
tator, annihilation and contributions from intermediate charmonia are explored in 
Section 0J In Section El a phenomenological discussion of the B K'-yy Dalitz 
region, in which the OPE is valid, and the di-photon invariant mass distributions is 
given in the SM and with NP. We summarize in Section IHl 

2 The lowest order h S77 amplitude 

We treat the effects of heavy degrees of freedom at the weak scale with an effective 
Hamiltonian known from b —>■ sj decays, e.g. [15] 





which contains four-fermion operators 



SaYLcp){Ci3jf,Lba) 



02 = ( 



SaYLCa){Cf3'y^,Lb/3), 



(3) 



and the electromagnetic dipole operator 



e 



mbSaO^^^RbaF, 



(4) 



167r2 
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Figure 1: Leading order Feynman diagrams contributing to b ^ S77 decays. The 
IPI diagram illustrates the insertion of four-fermion operators (left plot), while the 
IPR diagrams (right two plots) represent the insertion of 0-j. The indices i,j = 1,2 
(with i 7^ j) denote both photons and their interchanged counterparts. 



Here, a and /3 are color indices, L,R={1 ^ 75)/2 are chiral projectors and 
denotes the QED field strength tensor. In Eq. we neglect the penguin opera- 
tors s'jfj^Lb'^qj^L/Rq and double Cabibbo-suppressed contributions proportional 
to VubV*g. We further neglect the strange quark mass throughout this paper. An- 
alytical formulas for the Wilson coefficients Cj(/i) can be seen in [15]. Using the 
equations of motion, one can show that there are not further independent, gauge 
invariant operators with two photons and with dimension equal or less than six. 
Hence, the Hamiltonian in Eq. (j2)) is the same for 6 ^ 57 and b 577 decays. 

The lowest order Feynman diagrams contributing to 6 ^ 577 decays are shown 
in Fig ^ Only the sum of the two IPR diagrams is gauge invariant [2]. The 
corresponding amplitude can be written as [5] 

^(&-^77) = --^V,,V:Mp') ■ [{NcC^ + C2)QlnM'' (5) 
+ C7QdW7^'']-utip)e^ih)e,ik2), 

where p{p') represents the momentum of the 6(s)-quark. We denote by ki,k2 the 
4-momenta and by e(A;i,2) the polarization vectors of the photons. The auxiliary 
tensors are written as 



{p — ^2)^ — fnl {p' + k 



2 



+ {ki,^<-^ k2,i^) , (7) 



where Nc = 3 denotes the number of colors, Qu = 2/3 and Qd = —1/3 are the 
up-type and down-type quark electric charges. Note that 1^2*'^ is proportional to 
the Adler- Rosenberg tensor [16]. We further suppress the renormalization scale 
dependence of the Wilson coefficients. To the order we are working, the Wilson 
coefficients are evaluated at leading log at /i ~ nib. We further use e°^^^ = —1. The 
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loop function Kg is defined as [17] 

"'^ dx 

^9 JO 

arctan 



f^iQ ,mg) = 7^ + — — In (1 - ZgX + ZgX ) 



4:-Zo 



- + - 

2 ^ 2o 



— 2z7r In 



wliere = /Tn?„ and = {ki + ^2)^. It assumes tlie following limits 



if Zq<A , 



otherwise, 



-2^ + 0{qVm'g) for « m^, 

\ + 0{ml/q^,rnl/q^ X logs) for q"^ ml. 



(9) 



3 B ^ ^^77 decays 

In this section we obtain explicit expressions for the B K'j'y'^ decay distributions. 
A general parametrization of the amplitude and decay rate is given in Section 13.11 
We perform an expansion in inverse scales of order rrih in Section 13.21 to render 
the matrix element of the IPR contribution local. This enables us to express the 
B — > K'j'j matrix element in terms oi B —y K form factors. The resulting spectrum 
is presented in Section 13.31 



3.1 The general amplitude and decay spectra 

We write the general B — * 7^77 amplitude as 

AiB ^ 7^77) = T,,e^{h)e''{h), (10) 

where 

T,. = A{g^,-^^^)+iBe,,^pk^,4 (11) 



+ 


D 


k k 

Pb ■ kipBvk2fj, + Pb ■ k2PBfikiy Pb ■ kipB ■ ^2 , / ki 

fcl ■ K2 


■ k2PBiiPBu 


+ 


iC 


h 


kl ■ /i^2C/^!^Q/3(^2 ~l~ ^1 )Pb ~^ (^2/iCi/a/37 + ^l;/C/ja/37)^l ^2^5 




+ 


iC 




kl ■ ^2^/11^0/3(^2 ~ ^1 ~^ {k2fi^uaf3-y ~ ^li/C/io/37) ^1 ^2^*5 




+ 


iD 


f 


Pb ■ kie^yapk2P% + Pb ■ k2ef_iua/3kiPB + {PBfj.^uaf3-r + PBu^fiaP"/)kik2p'^ 


+ 


iD 




Pb ■ kie^uai3k2PB - Pb ■ k2e^uai3kiPB + {PBi^eua^-r PBuef,ai3^)k^k^p]^ 



"'"Since we are not concerned about CP violation and to avoid clutter "we use the convention 
B = bq throughout this "work. 
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which is manifestly gauge invariant. Here, we denote by pb the 4-momentum of 
the 5-meson. The form factors A,B,D, C^, are functions of two kinematical 
variables. We choose later to use the photon energies Ei, E2 defined in the B-meson 
rest frame or the invariant mass q^. The amplitude squared is then given as 

\AiB ^ K^-f)\^ = 2\A\'' - 2Re{AD*) {mlh ■ k2 - 2h ■ psh ■ Pb) (12) 
+ {\D\^ + 2|Z}+n (m|fci ■ k2 - 2ki ■ pBk2 ■ Pb? 
+ 2(A;i ■ k2f\B - mlD- + k^ ■ pb{C+ - C") - ^2 ■ Pb{C^ + C-)\\ 

and the double differential rates (which includes a factor of 1/2 for identical photons) 

dV |^(fi^ir77)p 



dq^dEi 256m|7r3 



(13) 



The full photon invariant mass spectrum dV/dq^ can be obtained by integration 
with 

Er"^"^ = ^{<+^'-<^ V^^^^'O) ' (14) 



where A(a, b, c) = + 6^ + — 2ab — 2ac — 2bc. We also use 

dT \A{B ^ K-f-f)\'^ 



dEidE2 USmsTi^ 



(15) 



3.2 The OPE 



The matrix element of -B — K'j'-f decays induced by the dipole operator O7 is a 
non-local one, see Fig. ^ and difficult to estimate model-independently. Here our 
strategy is to go to a kinematical region where both the internal s- and 6-quarks are 
far off-shell, with virtualities of order m^. To see that there exists such a region of 
phase space, we investigate the propagators of the IPR diagrams in Eq. (jH)), which 
are given as (in the bottom rest frame) 

{Ql? = {P + hf = ml- 2m,E, , {Q\f = ~ [{p - k^f - ml] = 2m,Ei, (16) 
{Q2? = (P + kif = ml- 2m,E2 , {Q\f = - [{p - ^2)' - m^] = 2m^E2. (17) 

For photon energies neither too soft nor too hard with energies near their maximal 
value all Ql are hard. For example, if both photons have energies rrih/S, then 
Q1.2 — "^fe/v^ and Q\2= \/2/3. Note that here the light quark has large energy 
mf,/3 as well, and q^ = ml/ 3. 

As a next step, we integrate out the large scales associated with the intermediate 
quark propagators. We construct the resulting local operators out of a bottom heavy 
HQET quark hy and a strange collinear SCET quark x [13,14]. Here, v = PB/f^iB 
and a light-like vector n = px/E^, where Ek,Pk denotes the energy, 4-momentum 
of the kaon. Hence, we perform an OPE in Aqcd/Q where Q = {mi,, Ex, Q12, \/q^}- 
In the following we neglect mx and mB — m^f,. 
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For the lowest order matching onto the b — > 577 amphtude, the following gauge 
invariant operators with dimension 8 with a heavy and a coUinear quark and two 
photons are needed after using the equations of motion (Wilson lines are understood 
in the definition of the field x) '■ 

Qi = '^x<y,u<yapRh.FfF^\ q; = '^x<y,u<yapRKFfFr , (18) 

Q2 = -2im^x(y^,RKF!^''Ffv''nP, Q!^ = -2imkX(y,..RhvF^'' FfV^'n^ , (19) 
Q3 = xYLKFfD^F^f,^ Q'^, = xi^LKFfD^F^^^, (20) 

where F^y = l/2e^i,Q,^F"^. The primed operators are related to the un-primed ones 
by interchanging photons. Then we obtain 

The matching of the IPI contribution from 4-Fermi operators Oi^2 is readily com- 
puted: 

sWfbeihfeihy ^ 2 ^^\^\ (22) 

Q 

Note that all IPR and IPI diagrams are of the same order in l/Q. This is in contrast 
to the behavior in S ^ 77 decays, where the IPR contribution with intermediate 
s-quark gives the leading power result [6]. 




3.3 The B K^^ matrix element and decay distribution 

The matrix elements of the local operators contributing to S ^ ^^77 decays are 
obtained from tree level matching of the QCD onto the SCET currents [13,18,19] as 

{K{n)\xK\B{v)) = 2EKaEK), (23) 

{K{n)\xiMB{v)) = 2EKC{EK)n^, (24) 

{K{n)\xa^MB{v)) = -2iEKCiEK) {v^n, - v^n^) , (25) 

{K{n)\x<J^,,-f5K\B{v)) = -2EKC{EK)e^,.apv''n'^, (26) 

where the form factor ( can be identified with the QCD form factor in the usual 
parametrization, see e.g. [20], as C = /+• This is the only form factor remaining in 
the symmetry limit, which enters all B ^ K transitions such as S — > K£'^£~ or 
B — > Kuu decays. 
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We obtain then the following B — * K^y^y amplitude in terms of the parametriza- 
tion given in Eq. (fTTj) as 

A = -KC{EK)QdmbCj-^{2Ei + 2E2 - mB){^E^E2 - E^m^ - E^mji), (27) 

B = -k2C{Ek) ^Q'i'^bCr ^sE^^E^^l + ^^'^^ + C2)QlK}j , (28) 

D = KQ{EK)QdmbC-j -— ^ , 29 

2EiE2{'mB — 2Eij[mB — 2E2)m^ 



Ei-E, 



C+ = tiaEK)Qdm,Cj^^-^^, (30) 

^ =-'^2^^^-)(2i^, + 2E,-m,)m,' ^2^) 

- .t(F \nmr iE,-E2){2E, + 2E2-ms) 
D - <iEK)Qdm,C, ^^^^^^^^^^ _ 2^^^^^^ _ ^^^y (32) 

D = 0, (33) 
where 

K = +lA-^VtbV*aem. (34) 
V ZTT 

Here we use rrib ^ rriB but keep explicitly a factor of mj, next to C7 from the definition 
of the dipole operator given in Eq. 

From Eqs. (jHI) and (fT^ it follows that the B K'j'j rate vanishes for ^ 0. 
This must be so since in the limit that one of the photon energies vanishes the double 
photon rate does because the decay B — > K'j is forbidden. As a further consequence, 
phase space integration of the B K'j'j amplitude squared is IR finite, as opposed 
to the inclusive B ^^77 decays, which need the cancellation with the virtual 
electromagnetic corrections in the low photon energy limit [5]. 

We obtain for the differential decay spectrum 

dT al^Gl\V,,V:f ^|2 r 2^2,^ , {ms~2{E, + E2)? 



X y24:ElEl - 47715^1^2(^1 + E2) + ml{El + E^ 

+ 'i2m,QdQlRe{CjKl){C^Nc + ^2)777^(777^ -E^- E2){mB - 2{E^ + E2)) 

+ ?>2QI\k,\\CiNc + C2)'777|(777B - E^ - ^2)^. (35) 



We recall that Eq. (j35|) is valid in the domain of the Dalitz plot with a fast light 



meson with energy Ek ^ ^qcd and Q'l'\ q'^ hard 
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4 Additional contributions 



In this section we discuss contributions to B ^ K'^y decays beyond Eq. (j35|l . These 
are if a photon gets radiated from the spectator quark, which is discussed in Section 
14.11 and annihilation contributions, see Section 14.21 Effects of this type can be 
experimentally accessed by comparing neutral versus charged B decays. In Section 
14.31 we further discuss contributions from intermediate hadronic states. 

4.1 Photon radiation off the spectator 

The OPE is invalidated by contributions where one of the photons is emitted by 
the spectator, for an example see the left diagram of Fig. |21 In the case of soft 
gluon exchanges the momenta of the final state s- and anti-quark are ps ~ —ki 
and pq ~ — ^2 up to order Aqcd in the decaying heavy meson rest frame. This 
configuration does not allow a Kaon to be formed if the angle between the photons 
is large sin 6' ^ \k±\/EK ^ Kqcd/'^b- This condition is indeed satisfied in the OPE 
region around Ei^2,K ~ ^b/^i where the 5-decay products form a Mercedes Benz 
star in the 5-rest frame, i.e. have angles of 120° between each other. Therefore, 
the exchange of soft gluons in these types of diagrams is excluded by the special 
kinematics. 

There are further hard scattering corrections from energetic gluon exchange, 
which are, however, suppressed by the strong coupling constant. They can produce 
the K also in a highly asymmetric configuration where Ps — Pk and are not power 
suppressed. We assume here that the endpoint suppression is sufficient to render 
this contribution subdominant. Note that due to the existence of a hard line in these 
diagrams, in the effective theory they get contributions from SCET operators [21] 
involving collinear gluons. For example, in the left diagram shown in Fig. |21 the 
s-quark emitted from the FCNC vertex has virtuality (Q^)^ ~ ^i^b) ^"^^ after 
integrating out these scales it is mediated by an operator with field content x^v plus 
one photon and a collinear gluon. Analogous statements hold for the other gluon 
exchange topologies which are not shown in Fig. |21 




Figure 2: Contributions to the B K'yy amplitude beyond the form factor-type 
OPE. Left hand side: radiation off the spectator, the blob denotes an insertion of 
the dipole operator O-j. Right hand side: Leading power annihilation contributions, 
the blob denotes a Four-Fermi operator and the second photon can be attached to 
any of the crosses. Diagrams with interchanged photons are not shown. 
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4.2 Annihilation topologies 



If we take additional operators in the effective Hamiltonian Eq. Q into account, 
annihilations topologies in B ^ decays arise. Weak annihilation (WA) con- 
tributions in charged -B^ K'^'-ff decays are induced by the operators 0'^2 — 
{s'y'^Lu){'LL'~f^Lb), which are obtained from Oi 2 by replacing c with u, and are CKM 
suppressed. In both neutral and charged B decays penguin annihilation is possible. 
It is suppressed by the small Wilson coefficients of the penguin operators which are 
of order IQ-^ in the SM. 

The leading annihilation diagrams are the ones with at least one photon emitted 
from the light anti-quark in the S-meson, see Fig. El Then one propagator is off- 
shell by order ~ QAqcd- Dimensional analysis indicates that this contribution is 
leading power with respect to the form factor-type contribution^ 

•^^•^^ oc 0(1) with fB oc m^^^'^, /+ oc m^^^"^, (36) 



mB^QCDf^ 



therefore it is only suppressed by small CKM elements or small Wilson coefficients. 
Such leading power annihilation contributions have also been found in i? — > Ki~^i~ 
decays [22]. However, their impact is bigger in the two-photon decays than in the 
semileptonic decays since the dominant Wilson coefficients associated with the form 
factor type contribution is larger in the latter. We estimate for B — > K'yy decays 
^WA^^f arm factor _ VubV*J (VtbV*) C2/C7, a correctiou of order ten percent at the 
amplitude level. 



4.3 Long-distance effects 

Resonance effects through B {rj^ 77)^^ where rj^j. = rj^rj'^rjc and B [K* — >■ 
K'j)'j have been discussed in [10,11]. In particular, the latter decay chain has 
been used to model the O^-type contribution in these works. Relevant here to the 
kinematical OPE window are only the charmonium resonances, since invariant di- 
photon masses of m^, m'^, are too low and the decays with virtual K* has one photon 
with near maximal energy ~ mB/2, which is excluded. To be specific, the rjc and the 
Xco,c2 have masses in the signal region and sufficient branching ratios given as [23-25] 

B{B ^ Kr],) X B{r], 77) ~ 4 • 10"^ (37) 
B{B ^ Kxco) X Bixco 77) ^ 2 • IQ-^ (38) 
B{B ^ Kxc2) X B{xc2 ^ 77) < 1 ■ 10-^ (39) 

We neglect in our study the Xc2 since the production through B Kxc2 has not 
been seen yet and its impact on i? ^ ^^^77 is at least an order of magnitude below 
the other two charmonia. 

^More precisely, the scaling law of the heavy-to- light form factor is /+ oc TOj, ^^^/(l — g^/m^)^ 
[18], which simplifies to the "standard" form at large recoil also at — 0{ml) as long as is 
order one away from . 
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To estimate how severe the long-distance contamination is, we assume Breit- 
Wigner form for the scalar and pseudoscalar X = rjc, Xco resonance amplitudes 



A{B K{X 77)) = A{B ^ XK)- ^— ^,A{X ^ 77), (40) 

— nix + ^"^xf {X) 

where 

AiVc ^ 77) = ia{r,,)F^,F'^'' , A{Xco ^ 77) = a{xco)F^uF^'' , (41) 
and take |a(X)| and \A{B —>■ XK)\ from the measured rates 

|^(X^77)P< ^^^^_\AiB^XK)\'\pK 



TiX) = ^ '-^ ^ , T(B KX) = ^ ' ■ (42) 



This phenomenological parametrization has problems because the information about 
the strong phase is lost, a residual g^-dependence in the partial amplitudes is ne- 
glected and there is double counting to some extend. However, the typical size of 
the effect is reproduced correctly. We use this model to see whether short-distance 
physics can be extracted from B K'yy decays in the OPE region. This is worked 
out numerically in Sectional 

The pollution from cc is also known from rare semileptonic decays, which get 
a sizable background through B (v|/(''"v) — > £+£~)i^(*) that has to be removed 
by kinematical cuts. If duality would be perfectly at work, the IPI diagram with 
the charm loop would include the contributions from the charmonium resonances. 
For ~ ml ^ OPE methods have been suggested for B K*£^i^ decays 
for a model- independent calculation of the rate [26]. For B —>■ K'-ff decays, this 
amounts at lowest order to expanding n[q^,m1) = K(g^,0) + 0{ml/q'^), that is, 
using He = 1/2, see Eq. Q. However, in the double photon decays we cannot have 
q^ large enough and hence are too close to the charm threshold to use this expansion. 



5 Phenomenology of B ^ K^^ decays 

We discuss the OPE allowed phase space and give estimates for the B K'yy 
branching ratio in the SM in Section I^TD and in Section 1^7^ with NP. 

5.1 Dalitz plot and SM predictions 

The B — s> ^^77 Dalitz plot in the Ei — E2 plane is shown in Fig. |S1 The surviving 
phase space after cuts is the shaded triangle. We use £'1,2 < 2 GeV (thin solid 
lines) which imply Ek > 1.3 GeV, and q"^ > (long-dashed line). As can be 

seen, other cuts required by our OPE analysis, such as a minimal angle between 
the photon momenta for example 6 > 45° (dotted line) and a low energy cut on the 
photon energies of 1 GeV (short-dashed lines), are ineffective in the presence of the 
previous ones. The reference point £^1,2, Ek = m^/S is in the surviving region. The 
corresponding g^-range is numerically 5.6 GeV^ < q^ < 14.6 GeV^. Note that the 
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E2 GeV 




Figure 3: B ^ K^y^y phase space in the plane of the photon energies Ei and E2 in 
GeV. The thick solid line is the phase space boundary, the thin solid, dotted, short- 
dashed, long-dashed lines denote the cuts Ei^2 < 2 GeV, 9 > 45°, 1 GeV < Ei^2 and 
> Also shown is the point Ei = E2 = ms/'i. 



exclusion of the low di-photon mass events is not a big loss in the event rate since 
the B Kjj is suppressed by small g^. 

For our numerical analysis, we evaluate the SM Wilson coefficients at leading 
log at /i = rrih and take the B K form factor ( = f+ from Light cone-QCD sum 
rules /+(0) = 0.319 [20], which is in good agreement with the corresponding lattice- 
QCD resuh [27]. To be specific, we use Ci = -0.25, C2 = 1.11 and Cj = -0.35. We 
integrate Eq. ^ between maa;(Er", 1.2 GeV) and mm(Ep^, 2 GeV), see Eq. dH)), 
to obtain the di-photon mass spectrum dT/dq"^, which is shown in Fig. |3] As can 
be seen, the SD part of the spectrum (solid curve) is completely hidden behind 
the ?7c and the Xco resonance contribution (dash-dotted curve). The integrated SD 
branching ratio in the SM with cuts is small, AB{B — > K'j'y)^^^ ~ 1 x 10~^. It 
has an uncertainty from the form factor of about 20 % and about 50 % from the 
renormalization scale dependence, when /x is varied between mb/2 and 2mb. The 
order of magnitude of the branching ratio is consistent with the findings of Ref. [28], 
where the IPI contribution has not been taken into account and with no cuts. 

We also show for comparison the contribution from IPI SD diagrams only (dashed 
curve), which gives a somewhat reduced branching ratio B{B K'y'y)g[f'^^^ ~ 
0.5 X 10~^. Even if we integrate this contribution over a larger phase space region, 
the resulting branching ratio is about two orders of magnitude smaller than the 
corresponding ones ~ fewxlO"^ given in [10] (where a very large value of the form 
factor is used) and [11]. Such large B K'j'j SD branching ratios are also in 
conffict with the one of the inclusive decays B{B — > ^^77) ~ 3.7 x 10~^ [5]. 

5.2 New physics 

To be comparable to the resonance contributions at least in some region of phase 
space the SD contribution has to be lifted by roughly one order of magnitude above 
the SM one. We investigate here whether there exists such type of NP that is at the 
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Figure 4: Di-photon mass distribution of B ^ i^77 decays with the cuts speci- 
fied in Section{^from the OPE. The solid (dashed) curve denote the SM pure SD 
contribution with (and without) the IPR terms. The dotted (double- dotted) curve 
corresponds to a NP scenario with Cp = —700 (+700) and = 0. The dash-dotted 
curve is the SM including the resonance contributions from the rjc and the Xco- 



same time not violating other data. 

Beyond the SM, contributions to the dipole operator Cy are tightly constrained 
by data on 6 — >• 57 decays. In particular, IC7I is fixed at the 30 percent level [29]. The 
QCD penguin operators have small Wilson coefficients in the SM C^m^ /C2 ~ 10~^. 
Even an enhancement of a factor of 100, which is problematic because of data on 
B Ktt and other rare decays would constitute only an 0{1) effect in 6 ^ S77. 
Hence, these types of NP are inefficient to enhance the SD contribution in B ^ K'j'y 
decays to be above the rjc, Xco background. 

Let us consider scalar/pseudoscalar operators with leptons i = e, 

O's^P) = ^sRbiij,)i . (43) 

They contribute to the matrix element of double radiative decays through a lepton 
loop, similar to the charm loop shown in Figure ^ The contribution w.r.t. the 
charm loop scales roughly as aem/(47r) x Cg^p-^/{Qf^C2). The Wilson coefficients of 
the muonic operators {i = fi) are constrained by B{Bs —>■ fJ^^fi~) data to be order 
one [29], hence way too small to raise the B K77 SD contribution above the 
resonance contributions. Assuming lepton universality, i.e., C^(p) oc m^, which is 
realized if the operators are induced for example by neutral Higgs exchange, the 
electron contribution is tiny and jC'LpJ < 30, which is still not large enough. If one 
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relaxes lepton universality, however, one obtains model-independently |C*5(p)| ~ 700 
from B{Bs —>■ r+r^) < 5% [30]. Sizable enhancements of the electron couplings are 
excluded by data on b ^ se+e^ decays. The operators Cjj-p) enter b s-y decays 
as a 2-loop electromagnetic radiative correction, hence give only tiny effects in the 
single photon channel. We further checked that the one-loop mixing of Og — Op 
onto Og, see [29] for details, which gives a correction to 6 — >• se'^e~ and b sfi~^fj,~ 
decays, is consistent with current data. Note that the decoupling of 6 — > S7 decays 
from NP in 6 — > S77 has also been discussed in the context of SUSY with broken 
R-parity in Bs —>■ '~f'~f and B —>■ Xg'yy decays in Ref. [31]. 

Explicit calculation of the r-loop^ gives the following corrections to Eqs. ()27j) - 
dnni) a,s A + 6Ar, B B + 6Br etc., where 

6A, = n^C{EK)C^sEKmr(^Kr{A + z,)-^y (44) 

SBr = K^aEK)C}.^il-2Kr), (45) 

and all other 6C^, SD^, 6D^ vanish. Note that the matrix elements of C^5(p) 
can be deduced from the one of the fierzed QCD penguin Og, see e.g. [4,5]. The 
impact of Cp = ±700, Cg = and all other Wilson coefficients assuming SM 
values is shown in Fig. HI (dotted and double-dotted curves). As can be seen, even 
maximal NP enhancement brings the SD distribution barely above the resonance 
background outside the rjc and Xco peaks. For Cp = —700, = we obtain 
AB{B _ft'77)^p^ ~ (1 — 2) X 10~^, an order of magnitude above the SM value. 



6 Conclusions 

We presented predictions for the exclusive B — > 7^77 decay. In particular, we 
calculated the SD contribution to the matrix element by means of an OPE in 1/Q, 
where Q is a combined scale of order rrib and contains the 6-quark mass, the energy 
of the final meson in the B rest frame, the photon invariant mass and the virtualities 
of the intermediate quarks in the IPR b —>■ 577 diagrams. Then, around the Dalitz 
region where all decay products have energy ~ rrib/S, the b — > S77 amplitude can 
be expressed at lowest order in terms of local dimension 8 operators made out of 
a heavy HQET and a collinear SCET field plus two photons. The resulting local 
matrix element is given in terms of a form factor known from B — ^ K£~^i~ decays. 

We found that the resulting SD branching ratio in the SM is small, order 10~^, 
which is at variance with previous estimates [10,11]. Such small branching ratios 
for double radiative decays are not a complete surprise. In fact, being of the same 
order in aem as semileptonic rare decays with branching ratios of fewxlO"^ [20], we 
expect 
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^For the OPE this requires the operators Q4 — mi,xRhvF ■ F and Q5 = nibxRhyF ■ F. 



14 



The double photon decays are substantially suppressed with respect to the semilep- 
tonic ones, since the di-lepton operators with large coefficients |Cf]^| ~ C(4 — 5) 
are not contributing to the former, see e.g. [15]. 

In our analysis we further took long-distance effects via B r]cK '-fyK and 
B — > XcqK — >■ 77-?^ into account and found that the SM SD contribution is not 
accessible behind the resonance peaks, see Fig. HI We also discussed additional con- 
tributions from photon radiation off the spectator and annihilation diagrams, which 
induce sizable theoretical uncertainties and exhibit increased complexity compared 
to the ones in B K*^ decays due to the second (energetic) photon. On the 
other hand, non-standard scalar/pseudoscalar couplings to taus can give enhanced 
B — >■ ^^77 branching ratios of the order 10~^. Note that such a NP scenario is con- 
sistent with data on other rare decays such as 6 — 57 and h si~^i~ . In particular, 
large couplings are allowed because FCNC decays into a tau pair are essentially un- 
constrained to date. In other words, decays such as B^ — > r+r" - and 5 ^ K^*^ 
have sizable room for NP. The corresponding maximal SD B 7^77 spectrum 
with NP in the r couplings leaks marginally outside the resonance background. We 
conclude that it is very difficult to test SD physics with exclusive B — > ^^77 decays. 

Very similar arguments apply to B —>■ K*^^ decays. We obtained for the pure 
SD SM estimate of the IPI induced contribution B{B —>■ K*'~f'~f) ~ few x 10^^, which 
is two orders of magnitude below the corresponding value quoted in [32]. 

Note added: By means of the Schouten identity [33] only four of the form 
factors in Eq. (fTT|) correspond to independent Lorentz structures. In fact, and 
D~ can be absorbed into B, such that redefining B ^ B — m\D^ + ki ■ Pb{C^ — 
C^) — k2 ■ Pb{C^ + C ) with the remaining form factors A, D, left unchanged 
gives identical results for the B ii'77 amplitude. 
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